On bipartite Rokhsar-Kivelson points and Cantor deconfinement 
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Quantum dimer models on bipartite lattices exhibit Rokhsar-Kivelson (RK) points with exactly 
known critical ground states and deconfined spinous. We examine generic, weak, perturbations 
around these points. In d = 2 + 1 we find a first order transition between a "plaquette" valence 
bond crystal and a region with a devil's staircase of commensurate and incommensurate valence bond 
crystals. In the part of the phase diagram where the staircase is incomplete, the incommensurate 
states exhibit a gapless photon and deconfined spinons on a set of finite measure, almost but not 
quite a deconfined phase in a compact U{1) gauge theory ind = 2-|-I! Ind = 3-|-lwe find 
a continuous transition between the U{\) resonating valence bond (RVB) phase and a deconfined 
staggered valence bond crystal. In an appendix we comment on analogous phenomena in quantum 
vertex models, most notably the existence of a continuous transition on the triangular lattice in 
d = 2 + l. 
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I. INTRODUCTION 

Quantum dimer models (QDMs) were introduced by 
Rokhsar and Kivelsoni to capture the low energy dynam- 
ics of valence bond dominated phases of quantum Heisen- 
berg antiferromagnets.^ Work on these models has estab- 
lished that their behavior on bipartite and non-bipartite 
lattices is fundamentally different. On the latter, they 
admit RVB phases with liquid dimer correlations and 
deconfined (test) spinons^*^ which have topological or- 
der of the Z2/BF varietyi^ On the former they admit 
only crystalline phasesiiSiSiiSiiiii^ (valence bond crys- 
tals) in d = 2, while supporting RVB phases of the U{1) 
or Coulomb variety in d > 2.^'^'^*-^^ 

One of the elegant features of the class of quantum 
dimer models introduced in Ref. Q is the existence of 
special points, christened RK points, at which the exact 
ground state wavefunctions are trivially determined and 
take the form of equal amplitude superpositions over sets 
of states connected by the dynamics — the classic RVB 
form. Additionally, test monomers have a vanishing in- 
teraction at these points, which translates into the de- 
confinement of gapped spinons in RVB languagei^ The 
appropriate phase diagram for the square lattice case is 
sketched in Fig. ^ 

In this paper we study quantum dimer models "near" 
the RK points on bipartite lattices, i.e. QDMs that con- 
sist of the RK point Hamiltonians plus small, but generic, 
perturbations. The restriction to bipartite lattices brings 
an important simplification. In a sense that we will make 
precise below, bipartite RK points are described by gap- 
less Gaussian field theories and this will greatly simplify 
the task of analyzing perturbations. 

Our basic new results are twofold. First, we establish 
the degree of criticality of the RK points — more precisely 
of the associated Gaussian field theories which are fixed 
points of the obvious renormalization group and to which 
we shall refer for clarity as the RK fixed points (RKFPs) . 
In d = 2 -|- 1 we show that on the honeycomb and square 



lattices the RKFPs describe multicritical points with two 
relevant symmetric operators. In d = 3 -|- 1, the cubic 
lattice RKFP describes a critical point with just one rel- 
evant symmetric operator. Second, we show that the 
purely "staggered" phase that borders the RK points in 
the simplest QDMs can be pushed a finite distance away 
in the perturbed models and for d = 2 -I- 1 the interme- 
diate region is filled in by a devil's staircase of commen- 
surate and incommensurate valence bond crystals which 
exhibit staggered order but also additional Bragg peaks. 
The commensurate crystals are gapped and confining. 
The incommensurate crystals exhibit a gapless photon 
(phason) and are deconfining, for sufficiently weak per- 
turbations. Further, although they occur at (boundary) 
points, they form a generalized Cantor set of finite mea- 
sure in the same limit. We refer to this phenomenon, 
which comes remarkably close to a deconfined phase in 
a compact U{1) gauge theory in d = 2 -|- 1, as Cantor 
deconfinement. In fact, near the RK point this region 
is practically a deconfined phase, as the commensurate 
confined phases occupy an asymptotically small fraction 
of the phase diagram in this limit. In addition, the gap 
in the confined phases is extremely small. The possibility 
of such phases was pointed out previously by Levito\iiS 
based on a mapping to the roughening problem. 

This paper is organized as follows. In Section ^ we 
discuss the height action for a 2 + 1-dimensional (gener- 
alized) quantum dimer models on the square and hon- 
eycomb lattices. In Section IIIII we discuss in detail the 
case of the honeycomb lattice and the mechanism that 
drives the quantum phase transition first order in that 
case. In Section IIVI we discuss the development of the 
tilted phase and in Section we show that due to the 
strong coupling nature of these generalized dimer mod- 
els, their tilted incommensurate phases generically ex- 
hibit a fully developed devil's staircase analogous to the 
one discussed in two-dimensional classical systems. In 
Section IVll we discuss the phase diagram for the case of 
the square lattice, and in Section IVIIII we discuss the 
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violation of the Landau rules in these quantum phase 
transitions, and the deconfining nature of these quantum 
critical points, including the relationship of RK points 
to the "deconfined critical points" discussed recently by 
Senthil and collaborators.^^ We summarize our conclu- 
sions in Section llXI In the appendix, we provide a brief 
discussion of analogous phenomena in quantum vertex 
models. Here, we find that the RKFP in c? = 2 + 1 that 
controls the transition between a crystalline state and a 
Cantor deconfined region is critical (rather than multi- 
critical) . 

II. HEIGHT ACTION IN d = 2 + 1 DIMENSIONS 

The Hamiltonian of the simplest quantum dimer model 
introduced in Ref. Q has the form: 

H = -t(p){\\\ + h.c)+v(p)(Z\ + \\\){\\\) 

H = ^t{\<j){€i\ + h.c) +v(]':>){':>\ + \a){a\){2.i) 

for the square and honeycomb lattices respectively. Here, 
a sum over all plaqucttes of the lattice is implicit. On 



the honeycomb latticejiSiii and with less confidence on 
the square lattice fi'i'SAiliiii it is known that the model 
exhibits three phases. 



For w ^ t it is in a "columnar" crystal phase which 
gives way by a first order transition to a "plaquette" 
phase as v/t is increased. The plaquette phase is char- 
acterized by an order parameter that vanishes contin- 
uously at w = t, the RK point. At the RK point there 
are multiple ground states — the exact ground state wave- 
functions are the equal amplitude superpositions of states 
in each sector of dimer configurations connected by the 
(off-diagonal) resonance term in Eq. H2.ll) . It is generally 
assumed that these sectors are defined by two winding 
numbers (see Ref. 15 for details) although a proof of er- 
godicity within sectors exists only for the zero winding 
sector>i2i For v > t the system enters a "staggered" phase 
where the ground states contain no flippable plaquettes. 
Cartoons of these three phases - columnar, plaquette and 
staggereed - for the square and honeycomb lattices can 
be found in Figs. 1 of Refs. Ull and lT2[ respectively. The 
square case is reproduced in Fig. ^ 




FIG. 1: Phase diagram of the square lattice quantum dimer model. In the top part of the figure, the mean tilt of the height 
surface is plotted as a function of v/t (see Section II). The corresponding dimer phases are sketched on the bottom part. The 
flat side (columnar and plaquette solids) terminate in the RK critical point {v/t — 1), which has no dimer long-range order. 
Here, for the RK quantum dimer model (dashed line), the tilt jumps discontinuously, corresponding to a first order transition 
into the staggerred solid. Upon inclusion of longer-ranged interactions, however, the tilt will ascend a devil's staircase from the 
fiat region, through a succession of incommensurate and commensurate phases. This phenomenon, (and the resulting structure 
factor) are discussed in detail in the text. 

I 



On bipartite lattices the hardcore dimer constraint can be converted into a Gauss law, V • E — 0. This is ac- 
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complished by endowing the links of the lattice with an 
orientation, so that they point from one sublattice to the 
other. Identifying an empty link with an electric flux — 1 
a link occupied by a dimer with a flux z — 1 (where z 
is the coordination of the lattice) one obtains a lattice 
electric field E with the desired property. In d = 2 + 1 
this can be solved by writing E = V x /i in terms of a 



scalar height field /t .— Here, the violation of the hard- 
core dimer covering constraint, through the presence of 
a monomer or overlapping dimers, shows up as a vortex 
in the height field. 

In terms of h, the properties of the above QDMs in 
the vicinity of their RK points are reproduced by the 
imaginary time Lagrangian 



^ - i(a,/i)2 + ip2(V/i)2 + ip4(v2/i)2 + Acos(2^/i) (2.2) 

with — p2 oc (v/t) + 1 changing sign precisely at the RK points. The last term keeps track of the discreteness of the 
microscopic heights, as is usual in height representations of problems in d = 2, and we will have more to say on that 
later. The square and honeycomb lattice problems differ in the relevant values of p4 which are (7r/32)^ (square lattice) 
and (tt/IS)^ (honeycomb). 

To determine the operator content of these theories, one also needs to specify the compactification radius of 
the height variables. This is determined by the minimal uniform shift of the heights which leads to an identical 
dimer configuration. For the case of the honeycomb and square lattices, the compactification radius equals 3 and 4, 
respectively. 

Finally, the details of the lattice structure enter in the relationship between lattice quantities and operators in the 
field theory. Specifically, long wavelength correlations of the dimer densities (where i labels the dimer direction) 
on the honeycomb lattice can be reconstructed by means of the leading operator identifications, 

ni — - ~ \dxh + i[exp(27ri/i/3) exp(47ria::/3) + h.c] 



1 1,1^, V3, ^, , 1, 
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712-- = -( — 5^, + — 9j,)/i+-[cxp(27ri/i/3)exp(47ria;/3 + 47ri/3) + h.c.] (2.3) 



^3 - i = i(-i5a; - ^a^)/! + i[exp(27ri/i/3)exp(47rzx/3 - 47ri/3) + h.c.]. 

In these expressions, x and y denote Cartesian coordinates, the x axis being perpendicular to one dimer direction. 
The unit of length is given by the separation of two neighbouring parallel dimers. The corresponding expressions on 
the square lattice are 



4 ~ 4^ ^ 2 
1 I, , 1 



-{-If+ydyh + -[(-!)" + h.c] (2.4) 



n^-- = -{-lT+y+^dxh+-[{-l)ycM^ih/2 + Tn/2) + h.c] 



The effective Lagrangian of Eq. H2.2|l was also consid- 
ered in Ref. |22 in the context of the theory of classical 
Lifshitz (multi)critical points in three dimensions and, 
for this reason, this effective theory has been dubbed 
the quantum Lifshitz model in Ref. This effective 
Lagrangian exhibits a three-dimensional analog of the 
physics of two-dimensional classical critical isotropic sys- 
tems. Thus, the action defined by Eq. H2.2|l exhibits a 
line of fixed points parametrized by pi with p2 = and 
A = 0, under the natural momentum and frequency shell 
renormalization group with z = 2. The RK points on 
both lattices flow under this renormalization group, to 
RKFPs on this line with the above values of /54. This ac- 
counts for their critical correlations and a mode spectrum 
uj ^ k'^. In the U{1) gauge theoretic interpretation of the 



dimer model, this mode is a photon with an anomalously 
soft dispersion. The irrelevance of the cosine is the nec- 
essary vanishing of the instanton effects that otherwise 
gap the spectrum of such compact confining theories24 
and this also accounts for the deconfined spinous. The 
height action further accounts for two other non-trivial 
features of the RK points, namely that they have degen- 
erate ground states in all winding number sectors (the 
RKFP point action is insensitive to gradients of heights) 
and that for the zero- winding states the equal time height 
correlations are logarithmic, 

(/i(x, r)/i(0, r)) = ln(^|x|/a) (2.5) 

and precisely those of the classical dimer problem. This 
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in turn allows us to fix the values of p4 quoted earlier 
as the classical dimer problem is exaclty solvable. (This 
also clarifies that these values are really the fixed point 
values.) 

On both lattices, the action fEa. l2.2|l accounts for the 
plaquette crystals for v < t a,s flat states of the height 
variable with different local fluctuations: when p2 > 0, A 
is relevant, i.e. it is dangerously irrelevant at the RKFPs. 
The selection of plaquette states indicates that A > for 
v/t < 1 on both lattices. For v > t, p2 < 0: the height 
tilts, and in the absence of any restoring term in Ea. l2.2l it 
attains its maximum tilt consistent with the microscopic 
constraints. This tilted phase translates into staggered 
dimer correlations which break translational symmetry 
on the square lattice and only rotational symmetry on 
the honeycomb. 

The transition between the plaquette and staggered 
phases is somewhat unusual. As p2 goes smoothly 
through zero at the RK point, the plaquette order van- 
ishes continuously there. Further, as v , there 
are two divergent length scales: ^ ^ \v/t — and 
^'^^^ = 6, 1 on the square and honeycomb 
lattices, respectively. The latter arises from the danger- 
ously irrelevant cosine and is the length scale on which 
the plaquette order appears; the value of follows from 
the dimension of the irrelevant cosine at the fixed point. 
However, the staggered phase emerges at full strength 
immediately and the ground state energy density has a 
derivative discontinuity at the RK point which indicates 
a first order transition in the thermodynamic classifica- 
tion. 

Let us comment briefly on the origins of the effective 
action of Eq. H2.2|l . The height representation for dimers 
was used by several workers on this problem; the recog- 
nition that it can be extended to capture the physics of 
the RK point is due to Henley and a more complete ac- 
count of the properties of the RK point and the transition 
about it was given in Ref. ,2. Henley's basic observation 
was that the Hamiltonian Eq. 12.11 also governs the mas- 
ter equation for the probability distribution of the classi- 
cal dimer problems endowed with the simplest plaquette 
flip dynamics.^- In height variables the classical dimer 
problems are in the rough phase and so it is reasonable 
that the dynamics is captured at long wavelengths by the 
Langevin equation 



drh{x) 



SS[h] 



+ C(x,t) 



(2.6) 



where S[h] ~ / (Pxiyh)"^ is the coarse grained entropy 
functional of the classical dimer problem. For Gaussian 
white noise distributed C, the average over trajectories is 
weighted by the Lagrangian, 



C^idrh 



6S[h] 
Sh{x) 



(2.7) 



which, upon dropping a total time derivative, is the ac- 
tion governing the RK points. 



This construction shows that the RK point action has 
higher dimension operators than one might have guessed 
on symmetry grounds alone. Indeed, the action H2.2(l 
doesn't distinguish our two lattices with their different 
symmetries. As advertised, we are interested in this pa- 
per in weak, but otherwise generic perturbations of the 
Hamiltonians H2.1|l about their RK points. Such generic 
perturbations of the RK point Hamiltonians will give rise 
to all operators allowed by symmetry and so we shall in- 
clude them in the analysis, with the proviso that they 
enter with small couplings. A second implication of con- 
sidering perturbed models will be that the values of the 
couplings in (|2.2() , in particular the marginal coupling p4 
will also be somewhat different from their values quoted 
earlier — these differences will again be small. In the fol- 
lowing we will refer to this set of modiflcations as working 
near the RKFPs without further explanation. 

A flnal comment is in order on an unusual aspect of the 
RKFPs. These fixed points, and indeed the entire fixed 
line, exhibit an enlarged symmetry in which the height 
field is free to tilt with any magnitude in any direction — 
which is then spontaneously broken by the ground states. 
While this appears to be non-generic, such extra symme- 
tries are a feature of any purely Gaussian fixed point. For 
instance at and above d = 4, the Ising transition is con- 
trolled by a Gaussian fixed point that is invariant under 
constant shifts of the scalar field. Much as in that case we 
work near zero field (magnetization), in our problem we 
will work near zero tilt. In both cases this is the correct 
choice for a sector of the space of couplings — the sector 
in which the transition is continuous in mean field theory. 



III. EXPANDED ACTION ON HONEYCOMB 
LATTICE 

We turn now to the honeycomb lattice. For honey- 
comb dimers, the heights live on the triangular lattice, 
which naively indicates an action rotationally invariant 
to fourth order in gradients. However, the microscopic 
heights have a sublattice structure: they take values 0,1 
and 2 modulo 3 on the three sublattices of the triangu- 
lar lattice. Consequently, the symmetries of rotation by 
7r/3 and inversion also involve sending h ~h. It fol- 
lows then that the allowed non-irrelevant operators at 
the RKFP include the relevant cubic term, 



.1 



V3, 



^3 = 53 {dxh){-dxh Ydyh){-dxh + —dyh) (3.1) 



in addition to the marginal coupling, 
£4 = 54 [V/i • \7hf . 



(3.2) 



The growth of 33 in the infrared indicates that the criti- 
cality of the RK point does not survive its inclusion — the 
form of the interaction indicates that when 53 is nonzero 
the system has already entered the tilted phase, although 
with the tilt now locked to the lattice. This leads us to 
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conclude that the transition between the flat plaquette 
phase and the tilted phase goes truly first order. 

If we tune to zero, we can examine the stability of 
the multicritical RK point. In a rcnormalization group 
treatment, where we keep a cutoff in space alone, the 
lowest non-trivial order flows of the two interactions are 



dX 




'dt ^ 




dg4 _ 


9 2 


dt 


A 3/2^4 : 



(3.3) 



whence both interactions are irrelevant for 94 > 0.-^- As 
(74 is marginally irrelevant it produces logarithmic correc- 
tions and the flows also renormalize p4 thus parking the 
system elsewhere on the fixed line. With these observa- 
tions we see that the multicritical behavior represented 
by the RK point is stable on a surface of codimension 2, 
provided the long wavelength p4 < (7r/4)^. The latter 
condition holds for small perturbations of Eq. H2.2|) . 

IV. THE TILTED PHASE 

Once (74 is large enough, the tilt no longer jumps to 
its maximum value on leaving the flat phase. When g^ is 
nonzero, the tilt prefers to point along one of the dimer 
directions or in between two dimer directions depending 
on the sign of this coupling. A similar role is played by 
g4 on the square lattice (see below). Here the restriction 
to small tilts is crucial: the states for g^ > and 53 < 
are not related by symmetry; in particular, at large tilts 
there are only states with one sign of the tilt, i.e. with the 
dimers dominantly oriented along the three bond direc- 
tions. We also note that on traversing the multicritical 
point with 173 — 0, the same orientational selection is now 
effected by the fifth order term, 

(4.1) 

which is thus dangerously irrelevant at the multicritical 
fixed point. 

The properties of the tilted phase can be analyzed by 
expanding in small fluctuations about the weakly tilted 
state (henceforth we assume 53 < 0, so that the tilt is 
along the x-axis, and that {g^l is small) 

/i(r, r) = C • r + Sti{r, r), C = C i. (4.2) 

The corresponding Goldstone action is given (to 
quadratic order in fluctuations) by 

(4.3) 

Higher order terms (in 5h) are irrelevant due to non- 
divergence of fluctuations. Provided \p2\ < 53/^4 then 



vt ~ vi and a single correlation length can be defined, 
~ vi/yj^, SO that for sufficiently large momenta, 
q > the spectrum is quantum critical, while as g ^ 
the modes acquire a stiffness, their dispersion is linear 
with longitudinal and transverse velocities given by vi 
and Vf. The average tilt of the height variable trans- 
lates into new Bragg peaks appearing in the structure 
factor for the dimer densities (as deduced from Eqs. 2.3 
and 2.4). They are of two types: first is the commen- 
surate Bragg peak at the wavevector of the maximally 
"staggered" state located at the origin for the honey- 
comb case (and at (tt, tt) in the case of the square lat- 
tice) , the second type are the incommensurate peaks dis- 
placed from the characteristic wavevectors of the colum- 
nar/plaquette pattern (e.g. (47r/3,0) and (tt, 0) on hon- 
eycomb and square lattices, respectively) by an amount 
proportional to C. In the tilted regime the intensity of 
the commensurate peak varies as , while that of the in- 
commensurate peaks as ^ (on the square lattice as 
^ i6yp4 Generally, as a result of quantum critical fluc- 
tuations, the two peaks differ in intensity parametrically. 
Both their intensities are for p4 corresponding to the 
the purely short range Hamiltonian of Eq. 12.11 The exis- 
tence of gapless modes has important implications, one of 
them is that test monomers (spinous) interact through a 
potential that generally grows only logarithmically with 
distance, as opposed to the linearly growing confining po- 
tential encountered in commensurate valence bond solids. 



V. INCOMMENSURATE LOCKED CRYSTAL 
PHASES: THE DEVIL'S STAIRCASE 

The actual state of the system involves an interplay 
between the propensity of the system to establish a tilt 
discussed thus far and the discreteness of the microscopic 
heights. As we will quantify a little later, the operators 
encoding the discreteness are much more irrelevant (al- 
though dangerously so, see below) at the RKFP than 
those responsible for the selection of the tilt magnitude 
and direction. Therefore, near the RKFP the problem 
can be treated in the order in which we have discussed 
it. Away from this limit one can worry about more com- 
plex phase diagrams. 

In order to treat height discreteness properly, we need 
to generalize our treatment of the discreteness of the 
height in Eq. (|2.2|) . There we included a potential that 
favors discrete fiat configurations of the height field but 
in doing so we ignored the possibility of the height locking 
into tilted configurations. To identify the full set of lock- 
ing potentials, we note that the height and the lattice 
points together define a three dimensional lattice with 
points (/i,x). A general locking potential, in particular 
the one generated from the microscopic constraints under 
rcnormalization, can be expanded in the set of reciprocal 
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lattice vectors {G} of this three-dimensional lattice, 

Viock(/i, x) = ^ e*(G.'>+G.-x) (5^1) 

{G} 

and thus each G defines a locking potential. Nicely 
enough, the height-space lattice for the honeycomb dimer 
problem is the simple cubic lattice with the [111] direc- 
tion measuring height, so working out the set of locking 
potentials is particularly simple. 

As discussed in the previous section, in the tilted phase 
the mean-square fluctuations of the height are finite. As 
a consequence of this, whenever the average tilt of the 
height surface is commensurate with the lattice, the cor- 
responding terms in the locking potential (above) are 
asymptotically relevant. At such points the ground state 
is a commensurate valence bond crystal with a gapped 
spectrum. At incommensurate tilts the outcome depends 
upon the strength of Viock relative to the remaining quan- 
tum fluctuations. When the latter are strong, as is the 
case for small tilts, the ground state at most tilts is 
an incommensurate valence bond crystal with a Bragg 
peak at the incommensurate wavevector but with a gap- 
less (phason) spectrum. In the gauge theoretic interpre- 
tation of the dimer model, this excitation is a photon 
and correspondingly test monomers are deconfined — in 
that they experience the logarithmic interaction of free 
two-dimensional electrodynamics. If the locking poten- 
tial is increased and/or the quantum fluctuations de- 
creased, we expect a version of Aubry's "breaking of 
analyticity"— transition beyond which the incommensu- 
rate ground states are pinned, their low-lying excitations 
are localized though still gapless, and the incommensu- 
rate ground states occupy a set of measure zero in the 
phase diagram. 

The weak locking regime exists, parametrically, at 
small tilts near the RK point. For concreteness, con- 
sider approaching the multicritical surface by tuning 
with p2 set equal to zero. The closer we get to the tran- 
sition, the higher order the commensuration that leads 
to locking, with |G| ~ 1/C or larger. This has various 
effects: first the coefficients of the corresponding lock- 
ing potentials generally decrease as one goes to larger G; 
second, the operators become increasingly irrelevant at 
the Rokhsar-Kivelson fixed point; and third, the length 
scale increases at the crossover from the quantum criti- 
cal regime where the locking potential is irrelevant to the 
tilted regime where it becomes relevant. These effects 
combine together to make the renormalized strength of 
the locking potential, and thus the size of the gap in each 
commensurate phase, as well as the width of the param- 
eter range in the phase diagram occupied by the phase, 
all vanish at small tilts C exponentially in 1/C^: e.g., 
the gap 

A ^ C"'/(p4C^) ^ (5,2) 
As the correlation length ^ cx 1/C, this is equivalent to 



an ordering length that grows rather fantastically, 

^ r'^'/"' (5.3) 

where a, a' are positive constants. Consequently, even 
the commensurate phases are for all practical purposes 
gapless in this weak-tilt regime. Thus when the tilting 
transition is either continuous or very weakly first-order, 
there is a deconfined, weakly-ordered dimer-crystal phase 
near the RK point where the ordering is generically in- 
commensurate! By contrast, the existence of the strong 
locking regime will need to be determined by actual so- 
lution of the relevant models — the existence of the fully 
tilted states is not proof of strong locking, merely of the 
existence of a maximum tilt as a consequence of the pe- 
riodicity of the height variable. 

Returning to the infinite volume and zero temperature 
phase diagram, the commensurate and incommensurate 
states will be interleaved in the classic manner of the 
devil's staircase studied in a variety of commensurate- 
incommensurate problemssSi At weak locking, close to 
the RK point, the staircase is incomplete, i.e. the de- 
confined incommensurate points will form a set of finite 
measure as the control parameter, p2 or 53, is varied. In 
fact, the fraction of the phase diagram occupied by the 
incommensurate deconfined states approaches unity as 
the RK point is approached. This is the phenomenon we 
have termed Cantor deconfinement. Farther away from 
the RK point the staircase could become complete, with 
the incommensurate points being a set of zero measure, 
before the maximally tilted phase is reached. Needless to 
say, this scenario assumes that no first order transitions 
intervene. 

Finally, a brief comment on Bragg peaks. The actual 
details of the ordering — periodic in the commensurate 
states and quasiperiodic in the incommensurate states — 
will cause the second peak to fragment into multiple and 
even a dense set of higher-order Bragg peaks. 

VI. EXPANDED ACTION ON THE SQUARE 
LATTICE 

On the square lattice the generic case requires that we 
(a) allow p4 to vary, (b) include the quadratic and strictly 
marginal term, 

>c™ = [(W (a2/,)2] , (6.1) 

and (c) include the interactions 

^int - 54 [V/i • Vhf + 54 [{d^hf + {dyh)^] ' . (6.2) 

The impact of (a) and (b) together is that we are inter- 
ested in flows about a two dimensional surface of fixed 
points. Unlike on the honeycomb lattice, there is no cu- 
bic invariant so that the transition is continuous at the 
level of mean field theory for the height action. How- 
ever, there are two tree-level marginal couplings whose 
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fate needs to be decided by fluctuations. Let us focus for 
now on the flow around the fixed fine with p4 = 0. 
The lowest non-trivial order flows of interest are: 

dg4 9 ,2 ~ 1 ~2n 

-dt = -M^'^^' + ^'^' + l''^ 

Interestingly, the flows in the ((74,34) plane are attracted 
to the origin only along the positive (74 axis; generically 
the flows run away to the region where the action is un- 
stable at quartic order in gradients. From this we con- 
clude that it is highly likely that the transition is driven 
first order by fluctuations in the generic case. Now the 
RK point behavior requires a further fine tuning, 54 = 0, 
and hence it represents the behavior of a multicritical 
surface of codimension two. We find that the inclusion 
of a small p4 does not alter the topology of the flows and 
there remains only one trajectory flowing into the flxed 
point at the origin. 

For the runaway flows, the pattern of symmetry break- 
ing is indicated by the initial sign of and depending on 
that we find four states with the tilt either aligned or at 
angle tt/A to the lattice axes. All of them exhibit, neces- 
sarily, a modulation of the dimer density at wavevector 
(tt, tt). If 54 = 0, the sixth order term is now responsible 
for the orientational selection. 

The analysis of the tilted phase now requires identify- 
ing the reciprocal lattice vectors of the diamond lattice, 
with its [100] direction measuring the height. The qual- 
itative features of the devil's staircase are exactly as in 
the honeycomb case. Bragg peaks now appear at (tt, tt) 
and near (tt, 0) and (0,7r). 

VII. THE CUBIC LATTICE 

We turn finally to the case of the d = 3 cubic lattice. 
Now the constraint is solved by introducing a gauge field 
A, in terms of which the action is of the form,"^^ 

£ = i(5*A)2 + ip2(V X Af + ^p^iV x V x A)^ (7.1) 

with p2 — again signaling the RK point. In gauge- 
invariant notation, the quadratic terms become the 
squares of the electric and the corresponding magnetic 
field ind = 3 + l. 

The discreteness of the microscopic fields is irrelevant 
in a finite interval around the RK point and hence can 
be ignored in its proximity. The only non-irrelevant term 
consistent with cubic symmetry missing from Eq. 17.11 is 

P4 5](<9.(V X A),)' , (7.2) 



which is strictly marginal. Hence the critical behavior 
represented by the RK point is (marginally) stable. The 
additional coupling leads to rotationally non-invariant 
critical points. 

For p2 < 0, we need to include the (dangerously) irrel- 
evant terms that control the magnitude of the "magnetic 
field" , B = V X A. These are and Bi + B^ + and 
lead to staggered crystals (now with a modulation of the 
dimer density at (tt, 7r,7r)) aligned parallel or intermedi- 
ate to the lattice axes. The connected correlations of the 
dimer density exhibit distorted dipolar correlations. This 
spectrum exhibits a gapless photon with fluctuations of 
the fields about their ground state values, and deconfined 
spinous. Altogether, one obtains critical points separat- 
ing a liquid of dimers from a crystal. At larger tilts con- 
finement and further symmetry breaking can be expected 
to set in, but that is outside the validity of our analysis. 

We also note that other related problems, such as the 
two-dimer model on the bipartite diamond lattice studied 
by Hermele et alM in its incarnation as a pyrochlore 
Ising antiferromagnet with Ising dynamics, will mutatis 
mutandis exhibit the same phase diagram. 



VIII. LANDAU RULES AND ORDER 
PARAMETER THEORIES 

In Ref.H three of us commented that the RK point in 
the simplest bipartite QDMs in d = 2 + 1 sat at the tran- 
sition between two symmetry incompatible valence bond 
crystals and attributed the divergence of the correlation 
length at it to the deconfinement of the gauge/height field 
at the RK point. Recently Senthil et alm^ have consid- 
ered a similar scenario for the transition between a Neel 
state and a valence bond crystal, although now with dy- 
namical spinous, and christened the intervening critical 
point a deconfined critical point. 

At issue in these discussions are two distinct but closely 
related possibilities. First, the violation of the Landau 
rule that a phase transition between two symmetry in- 
compatible phases must generically be first order, except 
at multicritical points. Second, the formulation of the 
theory of the phase transition in terms of fields other 
than the order parameters of the proximate phases. At a 
technical level this happens when the ordering is driven 
by dangerously irrelevant operators. 

As we have discussed in this paper, the situation near 
the RK points is more complicated than was assumed 
in Ref. z: Instead of a single tilting transition between 
the plaquette phase and the staggered phase there is a 
devil's staircase of commensurate incommensurate tran- 
sitions. Nevertheless, sufficiently close to the RK point, 
any commensurate ordering is extremely weak and for 
practical purposes one can think of the RKFPs control- 
ling a continuous transition between a flat phase with 
plaquette order and a tilted phase with staggered order as 
well as incommensurate order with a continuously vary- 
ing wavevector. The weakly tilted phase however is (log- 
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arithmically) deconfined, unlike the fully tilted staggered 
phase. 

In the quantum dimcr models analyzed in this paper, 
the RKFPs control multicritical surfaces — in contradic- 
tion with the claim in Ref. 2. However this is not in- 
evitable in this class of problems. A quantum vertex 
model discussed in the appendix exhibits a critical sur- 
face controlled by RKFPs. In either case the underlying 
height action (|2.2II is not derived from the order parame- 
ters of the proximate phases. In this sense they lie outside 
the Landau paradigm when applied to quantum phase 
transitions. 

Finally, we note that in the d = 3 -I- 1 cubic QDM, 
the tilting transition remains continuous for weakly per- 
turbed QDMs and the tilted phase exhibits only stag- 
gered correlations. This is now a continuous transition 
inside a deconfined region where Landau theory predicts 
a continuous transition. That said, the critical theory 
is still not what one would naively deduce from an order 
parameter analysis of the ordered (staggered) phase — the 
correlations of the order parameter at criticality have a 
dipolar form. The complication is the need to properly 
treat the local constraint V • E = 0, which is accounted 
for in the action Ea. l7.ll bv the introduction of the vector 
potential A. It is also the case that the ordered phase 
breaks rotational symmetry due to a dangerously irrel- 
evant operator. Altogether while the three dimensional 
problem isn't quite as unusual as the two dimensional 
one, it also furnishes an instance where the critical the- 
ory requires a departure from the standard cookbook. 



IX. SUMMARY 

We have analyzed the phase diagram in the vicinity 
of the RK points in the bipartite QDMs. In d = 2 
we find a first order transition separating confining pla- 
quette phases from a devil's staircase of commensu- 
rate/incommensurate and confining/deconfining valence 
bond crystals on both the honeycomb and square lattices 
with the RK points sitting on a multicritical surface of 
codimension two. In d = 3 the RKFP controls a con- 
tinuous transition between an RVB phase and a decon- 
fined staggered valence bond crystal. QDMs that exhibit 
these phenomena can be constructed by adding multi- 
dimer potential energies to the RK point Hamiltonians. 
Whether the phenomena here can be identified in a spin 
model or in other physical systems, is an interesting topic 
for future work. Finally, we note that as we were writing 
up this work there appeared Ref. |2^ which has consid- 
erable overlap with our work but which differs from our 
conclusions on some points. 
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APPENDIX A: QUANTUM VERTEX MODELS 

Quantum vertex models are close cousins of quantum 
dimer models. Here we begin with a Hilbert space la- 
belled by configurations of a classical vertex model, i.e. 
each bond hosts an Ising variable which we picture as an 
arrow marked on it. Then we introduce a local quantum 
dynamics which consists of reversing local closed loops of 
arrows. In cartoon form, the resonance move analogous 
to the one captured in Eg. 12. H is given by 



(Al) 



As the bond variables here are oriented their physical 
origins will neccessarily be distinct from those of dimers. 
For example, a quantum six vertex model was introduced 
in Ref. as a model for quantum effects in two di- 
mensional ice as well as via a mapping from a planar 
pyrchlore Ising magnet placed in a transverse magnetic 
field. Subsequently it also been argued to arise as an 
effective theory of the quantum fluctuations in the d- 
density wave state where the microscopic variables are 
oriented currents. "^^ 

An exhaustive analysis of the RK point manifold in 
the quantum eight vertex model is contained in Ref. IH- 
When critical, the model is equivalent to the six ver- 
tex model. In this case the conservation law at the ver- 
tices leads straightforwardly to a height representation in 
which the two sublattices of the square lattice host even 
and odd values of a height field with circumference two. 
A nice feature of the quantum six vertex problem is that 
one can find a family of RK points at which the value of 
/94 can be tuned continuously; the ground state wavefunc- 
tions in these cases are no longer equal amplitude super- 
positions but still exhibit isotropic critical correlations. 
The analysis of perturbations of these RK points now 
closely parallels our discussion of the square lattice QDM 
and we conclude that they govern a multicritical surface, 
likely between a plaquette phased and a devil's stair- 
case region. The significant differences from the dimcr 
case involve a) the degeneracies of the phases and b) 
the varying dimension of the vertex operator piece of the 
dimer operator. The latter causes the relative strengths 
of the staggered order and incommensurate order to vary 
in the tilted phase. For example, near the equal ampli- 
tude RK point, the staggered order is quadratic in the 
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tilt as against cubic for the incommensurate piece. 

Of the 32 possible vertices on a triangular lattice, 20 
respect the 'ice rule' of having the same number of ar- 
rows pointing in an out. The corresponding 20- vertex 
model (the ice model on the triangular lattice) thus per- 
mits a height representation. The heights live on the 
dual honeycomb lattice and are even and odd on its two 
sublattices with compactification radius two. The sym- 
metries now force the height action to be even in h and 
isotropic to fourth order in gradients. Consequently the 
line of RK fixed points is now stable with a dangerously 
irrelevant sixth order term picking one of six lattice di- 
rections to orient the tilted phase. The {h, x) lattice is 
an hep structure with in-plane bonds deleted and its re- 



ciprocal lattice vectors govern the commensurate states 
in the devil's staircase region. 

The values of p4 stemming from particular microscopic 
models in this family are under investigation and will 
be reported more fully elsewhere;'^^ here we content our- 
selves with noting that at the equal amplitude RK point 
the height is rough''^ and hence in its neighborhood we 
now have an example of critical behavior governed by RK 
fixed points and a Cantor deconfined region which can 
be accessed without dialing the extra parameter needed 
in the dimer models discussed previously in this paper. 
Based on experience with this class of models, it seems 
likely that the flat phase is again a plaquette phase which 
is symmetry incompatible with the Cantor region. 
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